function testDisKmeans()

DisKmeans_lambda = 1.0;
ClusterNum       =      3;
load('DisKdata', 'feat_set',  'gIdx');

sz      =   size(feat_set);
ss   =   sum(feat_set, 1)/sz(1);
feat_set    =   feat_set - repmat(ss, sz(1), 1);

%idx      =      CLUS(feat_set, DisKmeans_lambda, ClusterNum, gIdx);

Gram    =    feat_set*feat_set';
%idx2        =       CLUSKER(Gram, DisKmeans_lambda, ClusterNum, gIdx);

idx2        =       CLUSF(feat_set, DisKmeans_lambda, ClusterNum, gIdx);

%x =[1:2:7; 2:2:8]' ;
x   =   feat_set;
[U,S,V] = svd(x);%svd(x, 'econ');
